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ABSTRACT 

We use recent observations of the evolution of the type la Supernova (SN la) rate 
with redshift (Dahlcn e t al. 2Q03). the d ependence of the SN fa rate on the colours 
of the parent galaxies l)Mannucci et al . 2005), and the enhan cement of the SN fa 
rate in radio- loud early-type galaxies ((Delia Valle et al.1 12005|) to derive on robust 
empirical grounds, the distribution of the delay time (DTD) between the formation of 
the progenitor star and its explosion as a SN. Our analysis finds: i) delay times as long 
as 3-4 Gyr, derived from observations of SNe la at high redshift, cannot reproduce the 
dependence of the SN la rate on the colors and on the radio-luminosity of the parent 
galaxies, as observed in the local Universe; ii) the comparison between observed SN 
rates and a grid of theoretical "single-population" DTDs shows that only a few of 
them are possibly consistent with observations. The most successful models are all 
predicting a peak of SN explosions soon after star formation and an extended tail in 
the DTD, and can reproduce the data but only at a modest statistical confidence level; 
iii) present data are best matched by a bimodal DTD, in which about 50% of type 
la SNe (dubbed "prompt" SN la) explode soon after their stellar birth, in a time of 
the order of 10* years, while the remaining 50% ( "tardy" SN la) have a much wider 
distribution, well described by an exponential function with a decay time of about 3 
Gyr. 

The presence in the DTD of both a strong peak at early times and a prolonged 
exponential tail, coupled with the well established bimodal distribution of the decay 
rate (Amis) and the systematic difference observed in the expansion velocities of the 
ejecta of SNe la in Ellipticals and Spirals, suggests the existence of two classes of 
progenitors. We discuss the cosmological implications of this result and make simple 
predictions, which are testable with future instrumentation. 
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1 INTRODUCTION 

Type la supernovae (SNe) SNe are very important objects 
in modern cosmology because they are bright sources that 
can be detected up to large distances and it appears that 
their intrinsic luminosities can be inferred directly from their 
light curves. Exploiting these properties, the study of SN 
la at high redshi fts has allowed the dis c overy of the cos- 
mic acceleration (iPerlmutter et al 1 ll998l: iRiess et alJll998l: 
IPerlmutter et aLj^ggT ! Even though these objects are com- 
monly believed to be associated with the explosion of a de- 
generate star as a white dwarf (e.g. Hillebrandt & Niemeyer, 
2000) the nature of the progenitors of this class of SNe is 
not firmly established, and several explosion patterns are 
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possible (see, e.g.. Branch et al., 1995; Yungelson, 2004) and 
each of these may dominate at different redshifts. As a con- 
sequence, the existence of systematics affecting SNe la at 
different redshifts cannot be ruled out (e.g. Kobayashi et 
al., 1998; Nomoto et al., 2003) and it is worth being further 
investigated for possible cosmological implications. 

Different expl osi on m odels iGreggio fc Re nzini 
Yungelso n fc Livid l2000t iMatteucci fc Recchil 



1982 



2001 



Belczvnski et alJ l2005t iGreggid l2005^ predict different 



delay times between the formation of the progenitor system 
and the SN explosion. Differences in the expected delay 
times are testable by the observations (Madau, Delia Valle 
& Panagia 1998, Sadat et al. 1998, Dahlen & Fransson 1999) 
so that constraining the Delay Time Distribution (DTD) 
will permit one to ascertain the nature of SN la progenitors 
by confirming or excluding some of these models. 
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2 NEW OBSERVATIONAL EVIDENCE 

Recently, three important observational results were estab- 
lished: 



(i) The evolution of the SN la rat e with red- 



shift 


is now measured up to z~1.6 JHardin ct al. 


200C 


: iPain et alJ 120021: IStroleer 


12003: iMadewick ct al.. 


2002 


: ICaDDellaro et alJ 12004 


Gal Yam & Maod 12004 


Mannucci et alJl2005l: iBarris & Tonrj2006'). The results at 



the highest redshifts, derived by the GOODS collaboration 
JPahlen et all2004lstrolger et al.l2o'o 4. 2005) show that the 
SN rate rises up to z~0.8, when the Universe was 6.5 Gyr 
old (see panel b of Fig. 1), and decreases afterward. How- 
ever, the detection of SNe at z ^ 1.5 is very challenging for 
the current instrumentation, resulting in a small SN sample 
(2 SNe in the Dahlen et al. highest redshift bin), and in a 
very uncertain rate. 

This behavior can be compared with the cosmic Star 
Formation History (SFH) which continues to rise up to 
z~ 2.5, i.e., at a time abo ut 4 Gyr earlier (^adau ct al. 
Il99^ : 'Giavalisco et al."2004'). These results have been inter- 
preted bv .Dahlen et al.. (.200ll and lStrolger et alj (|200^ as 
evidence of a very long delay time (~ 3 — 4 Gyr) between the 
formation of the stars in the binary system and the explosion 
of a SN la. 

(ii) Recently we hav e measured the SN rat e per unit mass 
in the local Universe jMannucci et al.ll2005^ ■ and we have 
found a very strong dependence of the rate on the (B— K) 
colour of the parent galaxies: blue galaxies (latest Hubble 
types) exhibit a SN la rate a factor of ~ 30 higher than that 
of red galaxies (early types). This result indicates that the 
delay time must have a wide distribution. In star forming 
galaxies, the delay time must be as short as the timescale of 
colour evolution (~ 0.5 Gyr), while the existence of super- 
novae in galaxies without any recent star formation argues 
that a t least some SN la ha ve long delay times. For this 
reason annucci et al.' ^2005"), following earlier suggestions 
(Dallaoorta 1973; Delia Vallc & Livio 1994; Panagia 2000), 
have proposed, on robust empirical grounds, the existence of 
two populations of progenitors, one related to the young stel- 
lar population, with rates proportional to the recent SFR, 
the other related to the old stars, with rates proportional to 
the total stellar mass. 

(iii) Delia Valle et al. (2005) demonstrated that early- 
type radio-loud galaxies show a strong enhancement, by 
a factor of about 4, of the SN la rate with respect to 
the radio-quiet sample, with a constant rate between the 
sam ples being rejected at a 99.9 6% confidence level (see 
also iDella Valle fc Panaeial ll2003h ). Both the radio activ- 
ity and the SN rate enhancement are interpreted in terms of 
episodes of star formation due to merging with small galax- 
ies. Since each e pisode of radio activity is e stimate d to lasts 
about 10* years iSrinand fc Gopal-Krishnai.l99& : ,.Wan et alJ 
I2OOOI) , under this hypothesis the evolutionary time for most 
SN la in radio-loud galaxies must also be around 100 million 
years. 

All the above issues constitute observational links be- 
tween the epochs of star formation and SN explosion, and, 
therefore, can be used to constrain the DTD over differ- 
ent timescales. The evolution of the SN la rate with cosmic 
time is sensitive to long timescales (a few Gyr). The depen- 



dence of the local rate with the parent galaxy colour samples 
timescales of the order of the colour evolution of the galax- 
ies, i.e., 0.5-1 Gyr. The relation between SN rate and radio 
power gives information on the timescales of the order of 10® 
years, corresponding to the radio activity lifetime. 



3 PREDICTING THE SN RATES 

The aim of this work is to identify, among the possible 
DTDs, which ones can simultaneously match all these three 
results. We will test a set of possible DTDs spanning a broad 
range of possibilities, using both simple analytical formula- 
tions and the results from theoretical works. For each DTD 
we compute the corresponding expected rates, and compare 
the results with the above mentioned observational results. 

We have assumed that all stars in the 3-8 M© mass 
range are possible progenitors of SNe la (e.g., N omoto 
1994), and have a probability r] iMadau et al.lll993) of ac- 
tually e xploding as such w ithin an Hubble time from for- 
mation. iMadau et al.l (Il99^ have estimated this "explosion 
efficiency" to be of the order o f 5-10% . Similar results have 
been obtained bv iDahlen et alJ J2004) and Grcggio (2005|. 
We would like to note here that assuming a different mass 
range for the mass of the primary would only change the in- 
ferred overall explosion efficiency by a factor that depends on 
the adopted mass limits and the initial mass function (IMF) 
(see section 17^ without affecting the conclusions about the 
nature of the progenitors themselves. 



3.1 The evolution with redshift 

The evolution of the cosmic SN rate with redshift was mod- 
eled by a simple convolution of the assumed DTD with the 
observe d cosmic SFH corrected for dust extinction, as de- 
rived bv iGiavalisco et al.l (|20o3)- This is based on a decade 
of observations at different wavelengths and with various 
techniques, starting from the pioneering work of Steidel and 
coworkers ( Steidel ct al. 1996J . 

The need for a large correction for dust extinction in- 
troduce an important uncertainty. The SFH derived directly 
from UV observations is similar in shape to the corrected 
one, but is about a factor of 3 lower. As a result, the use 
of this SFH would produce very similar results but would 
imply a number of SNe per unit mass in stars three times 
larger. 

We adopted the functional form of the SFH as a 
fu nction of the cosmic time t (expressed in Gyr) as given 
bv lStroleer et al.l tQ04) . i.e. : 

SFH{t) = a (t'-e-'^" + d e'*(*-*°)/'=) Mg/yr/Mpc^ 

where a=0.182, b=1.26, c=1.865, d=0.071 and io is the cur- 
rent age of the universe, 13.47 Gyr for the assumed cosmol- 
ogy (h,^M,^A) = (0.7,0.3,0.7). 

To properly take into account the time elapsed between 
star formation and the creation of white dwarfs that even- 
tually will explode, the DTDs have been convolved with a 
function of time that, for a given stellar generation, gives the 
number of white dwarfs accumulating as a result of stellar 
evolution. Adopting a Salpeter Initial Mass Function (IMF), 
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half of the white dwarfs from stars in the 3-8 Mq mass range 
form within the first 130 Myr. 

For the c omparison, wc have used the observ ed rates 
inlDahlen et al . (2004). Strolger. (.2003.) and Mannucci et alJ 
i2005l) . The SN rate per unit mass of this latter work was 
transformed to a SN la rate p er unit volume b y using the 
local 2dF luminosity function dCole et al.ll200lh . obtaining 
3.5 • lO"'^ SN la/yr/Mpc^ 

3.2 The dependence on the colour of the parent 
galaxy 

In iMannucci et alJ ll2005fl we have already demonstrated 
that the dependence of the rate on the colour of the parent 
galaxy can be reproduced by a simplified model, which com- 
prises two components, one proportional to the SFR (about 
40% of the core-collapse SN rate) and one proportional to 
the total stellar mass, (about 4.4- 10~* SNe per year per 
10^° M0 of stars). This is consistent with a DTD having two 
components: a prompt one, exploding with a time scale com- 
parable to that of the core-collapse (CC) SNe; and a second 
component characterized by a long delay, so that the corre- 
sponding rate depends only on the integrated history of the 
parent galaxy. This result was obtained by using solely ob- 
served properties of galaxies and SNe, and, therefore, does 
not depend on any model of galaxy formation. 

Here, we adopt a more refined and accurate approach, 
and use galaxy formation models to link their observed prop- 
erties with their past history. We reproduce the properties of 
the galaxies by calculating a large number of model galaxies 
(see Tabl e m with widely different properties. We used the 
iBruzuarfc Chariot i2003l) models with different SF histories 
(from single burst to a rate extended over one full Hubble 
time), metallicities from 2% to 250% solar, multiple bursts 
of star formation, and ages between 30 Myr and 13 Gyr. For 
each model galaxy we compute the present-day (at z=0) 
(B— K) colour and the present-day SN la rate, obtained by 
convolving the SFH of that galaxy with the assumed DTD. 

The dependence of the rate on the colours is obtained 
by subdividing the galaxies into the same bins of colours, as 
done for the galaxy sample monitored for SN discoveries as 
in Mannucci et al. (2005), and averaging the SN rate in each 
bin. 

To compute a meaningful cosmic average to be com- 
pared with the observations, we impose that the collection 
of galaxy models used reproduces the main integrated prop- 
erties of the observed galaxy population. This was done by 
assigning appropriate weights to each galaxy model to best 
match both the observed distribution of (B— K) colour of 
the galaxies in the local universe and the evolution of the 
cosmic star formation rate. Many solutions are possible, cor- 
responding to different mixes of galaxy models. The vari- 
ance in these predictions was therefore estimated by creat- 
ing many different sets of galaxy samples that, with difi'erent 
mixes of galaxy models, can still satisfy the requirements on 
the (B— K) colours and the cosmic SFH with different mixes 
of models and different weights. In all cases the peak-to- 
peak uncertainty range, represented as a shaded area in the 
panels (c) of Figures 1-5, is small enough to conclude that 
these predictions are robust with respect to details of the 
galaxy models. 

For each model, we have computed the value of reduced 



Table 1 . List of the main properties of the used galaxy models 



Model SF history Metallicity 



modi 


0.1 Gyr long Burst 


z© 


mod2 


exp. 


decay 


witn T — u.z ^yr 


40% Zq 


mods 


exp. 


decay 


witn T — u.z Oyr 


Zq 


mod4 


exp. 


decay 


+ n n O r^irT. 

wn^n T — v.z \jyv 


250% Zq 


mods 


exp. 


decay 


with T=5.0 Gyr 


2% Zq 


mod6 


exp. 


decay 


with T=5.0 Gyr 


20% Zq 


mod7 


exp. 


decay 


with T=5.0 Gyr 


40% Zq 


mods 


exp. 


decay 


with T=5.0 Gyr 


Zq 


mod9 


exp. 


decay 


with T=8.0 Gyr 


20% Zq 


modlO 


exp. 


decay 


with T=8.0 Gyr 


40% Zq 


modll 


exp. 


decay 


with T=8.0 Gyr 


Zq 


modl2 


exp. 


decay 


with T=8.0 Gyr 


250% Zq 


modl3 


50% 


mods 


-1- 50% modi 


Zq 


modl4 


90% 


mods 


-1- 10% modi 


Zq 


modl5 


97% 


mods 


-1- 3% modi 


Zq 


modl6 


99% 


mods 


-1- 1% modi 


Zq 


modl7 


50% 


mods 


-1- 50% modi after 3 Gyr 


Zq 


modlS 


90% 


mods 


-1- 10% modi after 3 Gyr 


Zq 


modl9 


97% 


mods 


+ 3% modi after S Gyr 


Zq 


mod20 


99% 


mods 


-1- 1% modi after S Gyr 


Zq 


mod21 


50% 


mods 


-1- 50% modi after 10 Gyr 


Zq 


mod22 


90% 


mods 


-1- 10% modi after 10 Gyr 


Zq 


mod23 


97% 


mods 


-1- 3% modi after 10 Gyr 


Zq 


mod24 


99% 


mods 


-1- 1% modi after 10 Gyr 


Zq 



X^, i.e. Xr ~ X^/D.O.F. by fully considering the Poisson na- 
ture of the errors, and comparing the observed and expected 
numbers of the events. In the figures we report the associ- 
ated probability P that the discrepancies are merely due 
to statistical fluctuations, i.e., the probability that the ex- 
pected rates are consistent with the observed values. The 
reported values of Xr and P do not take into account the 
"model" uncertainties, shown as shaded areas in panels (c) 
and (d), as they are of systematic rather then statistical 
nature. Nevertheless, except for the case discussed in sec- 
tion 2] all the models contributing to the shaded areas give 
agreements comparable to the reported central value. 

3.3 The dependence on the radio power of the 
parent galaxy 

iDeUa VaUe fc Panaeial (|20o3) and iDella Valle et al.1 
have demonstrated that the SN la rate is higher in radio- 
loud Ellipticals than it is in radio-quiet ones (about a factor 
of 4). These authors proposed that both the radio activity 
and the excess in SN la rate are due to merging episodes 
with small companion galaxies: the accretion of material 
both fuels the central radio source and produces star for- 
mation, which boosts the SN la production. In this scenario 
the observed excess implies that a significant component of 
the SN la population is characterized by delay times as short 
as those associated with the radio activity. 

The lifetime of the radio emission after an episode 
of star formation in not well known, but several works 
( Sri nand fc Gopal-Krishna 199a : .Wan et al. 2000) con- 
straint it to be between 0.5 and 1.5- 10* yr. After this period 
the radio activity is expected to decay quickly. By assuming 
that the radio activity and the new episode of star forma- 
tion are coeval, the excess of SN la is directly related to 
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the integral of the DTD during the radio lifetime. To cal- 
culate the SN rate associated with a given radio power, we 
assume that, after a merging episode producing new stars 
for a few percent of the total mass, an early type galaxy 
is radio-loud for 0.5 — 1.5 • 10* yr, radio-faint for a similar 
amount of time, and ra dio-quiet afterward (see the defini- 
tions in lDella Valle et a l. (2005)). As a result, the SN rate of 
the radio-loud and radio- faint galaxies is expected to be pro- 
portional to the average of the DTD during the first ~ 10* 
and ~ 5 • 10* yr, respectively, while the SN rate of the radio- 
quiet galaxies is proportional to this average over one Hub- 
ble time. In other words, the DTDs of SNe in the radio-loud 
and radio-faint galaxies must show an excess at short times, 
otherwise the rate distribution with radio power would be 
flat. 



THE PREDICTIONS OF THE 
SINGLE-POPULATION DTDS 

Fitting exclu si vely the re dshift evolutio n of S N la rates 
Dahlen et alJ ll2004h and IStrolger et al.l l^2004^ (see also 
Strolger et al.l (120051) ') derived indications for long delay 
times, i.e., a DTD peaked at about 3.4 Gyr with basically 
no stars exploding during the first two Gyr. This result is 
the direct consequence of the observation of a decreasing 
rate at z>l. We have tested if such a DTD can reproduce 
the other two observational evidences. In this case the only 
parameter that can be varied to simultaneously reproduce 
the observations is the explosion efficiency r). The results 
are shown in Figure 1, where panel (a) shows the used 
DTD, panel (b) the evolution of the SN rate with redshift, 
and panels (c) and (d) the dependence of the SN rate with 
the (B— K) colour and radio power of the parent galaxies, 
respectively. It is apparent that such a DTD, while account- 
ing very well for the SN la rate evolution with redshift, 
is totally inadequate to explain the two trends observed 
in the local universe: the statistical probability that the 
expectations are in agreement with the observations are 
below 0.1% for the rate-colour plot and below 2.5% for the 
rate-radio power plot. This is true also when considering the 
uncertainties in the model, shown as a grey area in panel 
(c): in all cases the probability of agreement remains well 
below 0.1%. This demonstrates that the estimated delay 
times (about 3-4 Gyr) are far too long to reproduce the 
observed variation with radio power and the enhancement 
of the rate in the blue galaxies: as no SNe are expected to 
explode during the first couple of Gyr, this DTD cannot 
reproduce any dependence of the SN rate on the current 
SFR. 



Table 2. List of some of the tested models together with the 
resulting values of Xr- Iri boldface the values below 1.5. Many 
theoretical models providing poor matching of the data are not 
listed 



Looking for a DTD that satisfies all the constraints, 
we have investigated a large number of possible "single- 
population" models, i.e., DTDs that can be associated to 
a single progenitor population and be described by a sin- 
gle analytical law. We used DTDs characterized by different 
shapes (exponential decline, gaussian shape, and constant 
over one Hubble time) and characteristic times between 0.1 
and 6 Gyr. The results are summarized in the first part 
of Table 1. None of these simple DTDs can satisfy all of 
these observational constraints simultaneously. In all cases. 



Model 


V 




Xr 






(%) 


z 


colour 


radio 


Single population models 


Exponen. decay, t=1 Gyr 


5.3 


2.7 


7.2 


3.7 


Exponen. decay, t=2 Gyr 


4.0 


1.0 


4.5 


3.7 


Exponen. decay, t=3 Gyr 


3.3 


0.8 


4.5 


3.7 


Exponen. decay, t=6 Gyr 


3.2 


2.1 


15.7 


3.7 


Constant, 


2.5 


4.2 


29 


3.7 


Gauss at 0.05 Gyr, (t=0.01 Gyr 


6.5 


3.9 


32 


16 


Gauss at 0.5 Gyr, cr=0.1 Gyr 


6.0 


4.0 


37 


3.7 


Gauss at 1 Gyr, (7=0.2 Gyr 


5.4 


3.9 


25 


3.7 


Gauss at 1 Gyr, (7=1.0 Gyr 


4.6 


2.4 


9.2 


3.7 


Gauss at 2 Gyr, (7=0.4 Gyr 


4.0 


1.9 


23 


3.7 


Gauss at 2 Gyr, (7=2.0 Gyr 


3.7 


0.5 


15 


3.7 


Gauss at 3.4 Gyr, (7=0.68 Gyr 


3.6 


0.2 


22 


3.7 


Gauss at 4 Gyr, (7=2.0 Gyr 


3.1 


1.0 


14 


3.7 


Theoretical models 


YLOO" , DD-Ch, 


3.8 


1.0 


0.1 


2.4 


YLOO, He-ELD, 


4.5 


2.2 


34 


3.2 


YLOO, SG-Ch, 


5.0 


2.5 


21 


3.7 


YLOO, SG-ELD, 


4.0 


0.8 


11 


3.7 


BBR05^ SDS, aA=0.3 


4.0 


0.7 


1.2 


1.7 


BBR05, SWB, a\=1.0 


2.2 


3.5 


13 


3.7 


G05'= wide DD rnx=OA f3 = -0.9 


4.0 


1.1 


0.6 


3.1 


G05 close DD t„x=0.5 P = -0.75 


4.0 


1.0 


0.2 


2.9 


G05 SD Chandra 


4.1 


0.9 


0.6 


3.3 


G05 SD sub-chandra 


4.3 


1.2 


0.3 


3.0 


MROl'' SD 


4.7 


1.5 


1.6 


2.2 


Two Populations models 


50% prompt + 50% gauss. 4 Gyr 


4.3 


0.5 


5.0 


0.5 


50% prompt + 50% expon. 3 Gyr 


4.3 


1.2 


0.7 


0.4 


50% prompt + 50% const. 


3.0 


2.0 


6.2 


0.4 



"lYungelson fc Livid i200nh: 
i2005l) iMatteucci fc Recchi 



Belczvnski et al.l i200. 
20011) 



"iGreggiol 



at least one of the constraints is not satisfied and the devia- 
tions are much greater than the observational uncertainties. 

Within this class of models, the observations are best 
matched by an exponential distribution of delay times with 
e- folding time of 3 Gyr (see Figure 2) . This distribution pro- 
vides a rather good description of the observed rates as func- 
tion of redshift (panel b) and of the parent galaxy colours 
(panel c). However, this model is unable to describe satisfac- 
torily the variation of the rates with the radio-power of the 
parent galaxy (panel d). These results suggests that while a 
DTD that extends over several Gyr is needed, an additional 
contribution at early times (below 10* yr), is necessary to 
explain all observations. 



5 THEORETICAL MODELS 

Several authors have derived the DT D from the models 
of st e llar evolution and SN explosion j ^rgggi^^^R.engiJ 
19831: lYungelsori fc Liviol \2Q0(i IMatteucci fc R,ecchil l200lf 
" Belczvnski et alJl2005l: lGreggidl2005l) . In Table 2 we report 
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Figure 1. Properties of the SN rates derived for a DTD having 
a saussian shape, cen tered in 3.4 Gyr and with (7=0.68 Gyr, as in 
IStroIger et"ai] ^0^. Panel (a): the DTD itself, plotted as num- 
ber of explosions per unit time after star formation; Panel (b): 
the evolution of the rate along the cosmic age, with the upper 
axis showing the corresponding redshift. The solid line shows the 
result s for the consid ered DTD. Dat a are from [M aniiucci et al.i 
]2005[) (white circle) . Istrolgeil <200^ (starry dotT lDahlen et"] 
J2OO4) (black dots, with black error bars showing the Icr Poisson 
statistical errors and the gray bars the systematic uncertainties 
reported by the authors) . P is the statistical probability of agree- 
ment derived from the values and takes into account the sta- 
tistical errors only. Panel (c): the predictions (solid line and tri- 
angles) of this DTD for the dependence of the rates on the parent 
galaxy (B— K) colour, expressed in SNe per century per 10^" Mq 
of stellar mass, SNuM. The gray region shows the peak-to-peak 
uncertainties in these predictions, as explained in the text. The 
dots a nd the dashed line show the observations in lMannucci et all 
l|200fi) . Panel (d): variation of the rates in the early- type galaxies 
with the radio power of the parent galaxy. The black dots are the 
observed SN rates in iDella Valle et aTl l200fi) with 1-a Poisson 
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Figure 2. Same as Figure 1, but for a DTD that is an exponential 
function with an e- folding time of 3 Gyr. The rise of DTD at about 
10* yr in panel (a) is due to the building up of the white dwarf 
population. It is apparent that while the fit to the variation of 
the rates with redshift (panel b) host galaxy colours (panel c) 
are quite satisfactory when considering the uncertainties of the 
models, the behavior of the SN la rates with radio luminosity in 
elliptical galaxies (panel d) cannot be reproduced. 



the results of the comparison between "theory" and "obser- 
vations" by changing the only free parameter r]. All models 
predicting a narrow DTD are ruled out by observations. 

Some of the single-degenerate and double-degenerate 
models, which predict very broad DTDs, spanning delay 
times from 10^ yrs up to 10^" yrs. give interesting results. 
This is the case for a number ofjO reggio (2005) models^ both 
SD and DD (see e.g. Figure 3), the lYungelson fc Livid ll200d) 
DD Chandrasekar mass m odel (see Figure 4), and the the 
iMatteucci fc Recchil i200ll) SD model shown in Figure 5. In 
all these models the DTD peaks at about 0.6 — 2-10* yrs and 
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then decays rapidly, roughly like , becoming ten times 
lower after about 1 Gyr. The dependence of the rates with 
redshift and galaxy colours are satisfactory reproduced, even 
if in some cases the fast evolution tends to under-predict the 
SN rate in the reddest galaxies. Most of these DTDs predict 
that about 5-15% of the SNe explode within the first 10* yrs. 
As a consequence they produce a SN la rate in radio loud 
galaxies only 10-40% larger then in radio-quiet galaxies, in- 
stead of the observed factor of 4. Given the low number of 
observed events, we cannot exclude these models with a very 
high statistical confid ence level (Prad 8o=2-10%). As an illus- 
tration, let's consider iGreggiol l^20n!^^ DD model (Figure 3) 
that well reproduces panels (b) and (c) . In this case the peak 
of the DTD occurs at about 2 • 10* yrs and the fraction of 
SNe in the first 10* yrs is only 5%. As a consequence the 
enhancement of the rate in the radio-loud galaxies is only 
10%. 

Even if these models cannot be ruled out with an high 
degree of confidence, it is evident that a DTD with both 
more SNe at early time and a slower evolution afterward is 
needed to fully account for the observations. 



6 TWO POPULATIONS DTDS 

For these reasons we considered a set of "two populations" 
models in which the DTD is obtained as the sum of two 
distinct functions. In all cases, we added a "prompt" gaus- 
sian centered at 5 x lO'^ yr to a much slower function, either 
another gaussian or an exponentially declining function. We 
will refer to the former component as "prompt" exploders 
and to latter as "tardy" ones^ 

Figure 4 shows the results of a model in which 50% 
of the SNe derive from the "prompt" population, and the 
remaining 50% from the "tardy" one that consists in an 
exponentially declining function with characteristic time 
of 3 Gyr. These results shows that the observational data 
are better reproduced by a DTD with a peak at short 
times (below 10* yr) that includes about half of the SN 
la events, and an extension toward very long times, say, 
3 Gyr and beyond. Provided that the "tardy" component 
extends well beyond 3 Gyr, its shape is not well con- 
strained by the fit: exponential decays with characteristic 
times between 2.5 and 8 Gyr can still provide reasonable 
fits. These uncertainties will reduce considerably when the 
SN la rates at z>0.5 will be measured with greater accuracy. 

We notice that a bimodal distribution of delay times 
should not be regarded as a simple heuristic method to bet- 
ter fit the data, because there are theoretical models which 
actually predict a bimodal DTD. One of the best exam- 
ples, see Figure 7, is obtained by iBelczvnski et all (|200.4l 
with a single-degenerate model with reduced common enve- 
lope efficiency (aA = 0.3). This model predicts a bimodal 
DTD which peaks at 10* and 3 • 10® yr, and is due to the 
explosions of both He and C-O white dwarfs. This model 
correctly reproduces the evolution of the rate with redshift 



^ We chose to use "tardy", rather than "delayed" , to avoid any 
misunderstanding with the "delayed" detonation model adopted 
for type la supernovae, e.g. Woosley (1990), or more recently, 
Golombek & Niemeyer (2005) 
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Figure 3 . Same a s Figu re 1, but for the theoretical DTD pre- 
dicted by iGreggiJ {2003) for her DD model for a wide binary. 
The grey area in panel (d) corresponds to different assumptions 
for the lifetimes of the radio activity, from 0.5 to 1.5-10* yrs for the 
radio-loud and radio-faint phases, and to a fraction of new stars 
between 1 and 3% in mass. This DTD reproduces the observed 
trends in panels (b) and (c), while the expected enhancement of 
the rate in radio-loud galaxies is much lower than observed factor 
of 4. 



and its dependence on the colours, and it accounts for the 
enhancement in the radio-loud galaxies only qualitatively, 
because its "prompt" peak is centered at 10* yr instead at 
the best-fitting value of 5 ■ 10*^ yr. 

Bimodal DTDs can be produced also by models with 
more than one type of progenitors, for example in which 
both the single-degenerate and double-degenerate channels 
are active (see, for example, Nomoto et al., 2003). A bi- 
modal DTD is al s o nat urally produced by the SD model by 
iKobavashi et al] il998l) in which two different companion 
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Figure 4. Same as Figure 1, b ut for the DTD from the DD 
model with Chandrasekar mass by lYungelson fc Livid i2000D . 
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Figure 5. Same as Figur e 1, but for the theoretical DTD pre- 
dicted by the SD model in lMatteucci fc Recchil i200jl . 



stars are present: either a red giant with initial mass of 
about 1 Mq and orbital periods of tens to hundreds days, 
or a main-sequence star with mass ~ 2 — 3 Mq and periods 
of the order of a day. 



7 DISCUSSION 

7.1 Other evidences for bimodality 

We have shown that the observational constraints to the 
SN la rates, namely the rate evolution with redshift, the 
dependence of SN la rates with host galaxy colors, and the 
marked increase of SN la rates in radio-loud Ellipticals, are 
best reproduced if about half of the SNe explode within 10® 
yr from star formation ( "prompt" component) while the rest 



have explosion timescales of a few Gyr ("tardy" component). 
Particularly from our analysis a number of facts emerge: 

• After constraining the delay time distributions to match 
both the observations at high redshift and in the local uni- 
verse, we derive a DTD characterized by a strong peak at 
early times (~ 0.1 Mpc) followed by an exponential function 
with a decay time of about 3 Gyr. Our result rules out all 
models predicting narrow DTD, pure exponential or pure 
gaussian DTD and constant DTD. 

• There are some single-population models, characterized 
by broad distributions of delay times (3 orders of magni- 
tude), which can account the observations, but only at a 
modest statistical confidence level. Given the current uncer- 
tainties on the SN rates, they cannot be confidently ruled 
out. 

• The best match to the data is obtained with a bimodal 
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Figure 6. Same as Figure 1, but for a DTD constituted by 
equally contributions (50%) of an exponentially declining func- 
tion with characteristic time of 3 Gyr and a gaussian centered at 
5 X 10^ yr and ct = lO'^ yr. In panel (b) the dotted and dashed 
lines show, respectively, the contributions from the "prompt" and 
"tardy" components. The DTD simultaneously reproduces all the 
SN la rate observations. 



DTD. However, the existence of a bimodality in the DTD 
does not directly imply the existence of two distinct pop- 
ulations of SNe. This bimodality might be related to the 
bimodal distribution of Amis (the magnitude decline of the 
light curve between its maximum and 15 days later) exhib- 
ited by SNe la in Ellipticals and in spirals (Delia Valle et al. 
2005, their fig. 6, see also Altavilla et al. 2004). If this is the 
case, the bimodality of the DTD gives support to earlier sug- 
gestions which interpreted the bimodal behavior of Amis in 
terms of an age sequence (Ruiz-Lapuente et al. 1995, Hamuy 
et al. 1996, Howell 2001, van den Berg et al. 2005), with 
metallicity being less important (Ivanov et al. 2000, but see 
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Figure 7. Same as Figure 1, but for the theoretical DTD pre- 
dicted by Belczynski et al. (2005) for the single-degenerate model 
with reduced common envelope efficiency {a\ = 0.3). The best- 
fitting value of the total efficiency is t)=4.0%. In Panel (b) the 
dotted and dashed lines show, respectively, the contributions from 
the "prompt" and "tardy" components associated to the two 
peaks of the DTD. This bimodal DTD appears to reproduce rea- 
sonably well the SN la rate observations, although the position 
and the intensity of the "prompt" peak produce only a small en- 
hancement of the rate in the radio-loud Ellipticals. 



Kobayashi et al. 1998). In this scenario it is likely to expect 
that the bright SN events (with small Amis) are prefer- 
entially associated with the "prompt" population, which is 
more common in spiral and irregular galaxies (Hamuy et 
al. 2000), while the faint and rapidly decaying objects are 
more easily found in the "tardy" component, which is more 
common in early-type galaxies. If this is true, then we ex- 
pect that SN la in radio-loud Ellipticals, where the "prompt" 
population is expected to dominate, show Amis smaller that 
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those in the radio-quiet EUipticals. This issue is discussed in 
iDella Valle et alj J2005l) (see their fig. 6). The small number 
of events with published value of Amis prevents us from 
reaching any definite conclusion. 

• Finally we notice that the existence of a continuous 
sequence in the spectroscopic properties of maximum light 
spectra of SNe la (Branch et al. 2006) is not at variance with 
the existence of two different evolutionary paths to rise SNe 
la, as our results may suggest. What is relevant for our anal- 
ysis is how the spectroscopic properties of SNe la correlate 
with the stellar population from which the SN progenitors 
originate. The systematic differences measured in the expan- 
sion velocities of the ejecta of SNe la occurring in Ellipticals 
and Spirals (Branch & van den Bergh 1993; Benetti et al. 
2005) are indeed fully consistent with our results. 

7.2 Bimodality and stellar mass 

We note that the main sequence lifetime is about 5 x 10* 
yr for a star of SMq, 10* yr for 5.5Mq and about 4 x lO'^ 
yr for SMq (e.g., Girardi et al. 2000). Therefore, the SNe 
of the "prompt" peak, which include about 50% of the to- 
tal number of SN la events and explode within 10* yr from 
star formation, must all derive from stars with masses above 
5.5M0. Also, for a Salpeter IMF, the number of stars be- 
tween 5.5 and 8Mq are about a third of those between 3 
and 5.5M0. This implies that the SN la efficiency for higher 
mass progenitors (5.5-8M0) is about 3 times higher that 
for lower mass progenitors (3-5.5M0). Therefore, given an 
overall efficiency of 4.5% (see Figure 6), it follows that the ef- 
ficiency for higher mass stars is ?7(5.5-8M0)~6.8%, and the 
one for lower mass stars is 77(3-5. 5M0)~2. 3%. Therefore, 
the requirement that about 50% of the SN la must explode 
within 10* yrs is that the efficiency and the characteristic 
delay time are expected to change considerably at ^5.5Mq. 

It is important to realize that the explosion efficiency of 
the "prompt" SN la is determined unambiguously by their 
number and the mass range of the progenitors. On the other 
hand, if the remaining 50% SNIa, i.e. the "tardy" compo- 
nent, could arise not only from stars in the range 3-5. 5M0, 
but also from stars of lower masses, then their inferred ex- 
plosion efficiency for low mass stars would also decrease, 
because the available pool of stars would increase whereas 
the number of "tardy" SNe does not. 

We cannot draw conclusions on whether this change of 
the efficiency at about 5.5M0 is due to a different physical 
process (e.g., SD vs. DD) or to one and the same process op- 
erating in separate regions of the parameter space (e.g., sys- 
tematic differences of the binary systems as a function of the 
stellar mass). For example, the theoreti cal models produc- 
ing wide DTDs described in section |^ l|Yungelson fc Livid 
I2OOOI: iBelczvnski et alj l2005l : lGreggiolr2005ll can be modi- 
fied to reproduce the empirical evidence if this strong varia- 
tion of the efficiency with mass is introduced in the models, 
boosting the "prompt" part of the DTD. Thus, it could be 
that the binary fraction for primary stars with masses above 
5.5M0 is markedly higher than for lower mass stars. Or it 
could be that the distribution of secondary star masses is 
more skewed toward masses close to the primary star mass 
and, therefore, the mass transfer be more efficient and faster 
(see Pinsonneault & Stane k 2006 for a discussion). 

The existing models (iGreggid l2005l : IBelczvnski et all 



I2OO5I: iNomoto et alJliool) are not able to resolve this am- 
biguity because of both uncertainties in the model assump- 
tions and possible coexistence of different physical processes. 
However, we are confident that a judicious analysis of data 
obtained for a large sample of SN la over a suitably wide in- 
terval of redshifts will make it possible to clarify this issue. 

7.3 Consequences of the bimodality on Cosmology 

In addition to providing essential clues to the nature of SN 
la progenitors, our results have also important implications 
for cosmology: 

• The fraction of SNe coming from either of the two pop- 
ulations changes with cosmic time, as can be seen from Fig- 
ure 6: the "tardy" SNe dominate at z<1.3 and the "prompt" 
SNe above this limit. The ratio of the "prompt" SN rate to 
that of the "tardy" SNe changes from 0.5 in the local uni- 
verse to about 1.2 at z=1.5. Simi lar results are obtained 
for the single-degenerate model bv IBelczvnski et alj (l2005h 
(Figure 7). 

• It is reasonable to expect that the two populations of 
SNe can be distinguished also by some intrinsic properties. 
As an example, it is possible that "prompt" SNe la are, on 
average, more affected by dust extinction than the "tardy" 
compo nent, as they must explode closer to the formation 
cloud dSuUivan et al.ll2003ll . In this case, the average proper- 
ties of SNe la are expected to change with redshift, especially 
at z> 1 when the "prompt" SNe become common. The Hub- 
ble diagram s used to derive in formation on the cosmological 
parameters jRiess et alJl2004^ are, up to now, mostly based 
on SNe at z< 1 and, therefore, are expected to be domi- 
nated by the "tardy" population. As the ratio between the 
two different flavors of type SNe la changes with cosmic 
time, evolutionary effects should become important at high 
redshifts (Livio & Riess, 2006). 

• The l umino sit y-decline rate r elation for SN la 

JPskovskiil I1977I: jPhiUipsI I1993I: iHamuv et all Il996l: 
iPhiUips et alJ Il997^ is derived in the local universe 
and, therefore, under this scenario, is dominated by the 
"tardy" SNe. The evidence for a cosmological acceleration 
relies on the assumption that the same relation holds also 
at high redshift (see Rowan-Robinson 2002 and Leibundgut 
2004 for a discussion) . If the two populations follow slightly 
different relations, a bias is expected to emerge as a func- 
tion of redshift, especially when the "prompt" population 
becomes dominating, at z ~ 1.2. Thus, a confident use of 
SNe la for cosmology measurements at z> 1 would require 
a good understanding of the differences in properties of the 
two populations. 

7.4 Bimodality and metallicity evolution 

The existence of two populations of SN la has direct conse- 
quences also on the chemical evolution of the universe. 

• "Prompt" type la, having a redshift distribution simi- 
lar to the CC SNe, dominate the SN la population at high 
redshifts. Therefore, in the early Universe the production of 
Fe is expected to follow that of Oxygen, and the O/Fe abun- 
dance ratio should be relatively constant but appreciably 
higher than in the local Universe. When the SN la "tardy" 
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Figure 8. Ratio of the rates of the CC to la SNe as a func- 
tion of the redshift. The white and black do ts are observed values 
jMannucci et al.ll2005l: iDahlen et al.ll2004l). The lines show the 
predictions of the gaussian "sin gle-population" model in figure 1 
(dashed), the lYungelson fc Livig t200Q') model in figure 4 (dot- 
ted), and the "two-populations" in figure 6 (solid). The predic- 
tions use a Salpeter IMF and mass ranges of S-SMq (SNe la) and 
8-40Mq (CC SNe), and are scaled to match the observed values. 

component starts dominating, i.e. past the SFH peak at 
z~ 2, the Fe production is boosted and the 0/Fe ratio is 
expected to decrease rapidly to approach the "solar" values 
around redshifts <~ 0.5. Those aspects have recently been 
discussed in some detail by Scannapicco & Bildstcn ( 2005 ) 
who for their model calculations adopted t he simplified de- 
script ion of the SN la rates as derived bv iMannucci et aU 
(|2005|) in terms of a component proportional to the star for- 
mation rate (SFR) and another one that is proportional to 
the total stellar mass. 

• It is known that the intra-cluster medium is relatively 
rich in iron ([Fe/H]~ —0.5) and that the metallicity shows 
a very mild evolution with redshift (Tozzi et al., 2003). The 
observed iron mass is about a factor of 6 l arger than could 
have been produced by core-collapse SNe f|Maoz fc Gal-Yaml 
|2004) and a factor of 10 larger than that produced by the 
current rate of SN la (Renzini, 2004). The "two populations" 
model naturally explains these observations, as the current 
type la rate is just the long-time declining tail of a SN dis- 
tribution that peaked at early cosmic times. The amount 
of observed iron and its redshift evolution is reproduced by 
assuming an average age of the stars in clusters of 10 Gyr 
(see Matteucci et al., in preparation). 

7.5 Predictions 

The existence of the "prompt" and "tardy" populations of 
type la SNe can be tested by two observations: 

• The SN la rate is expected not to decrease significantly 
moving toward high redshifts up to z~2, at which the cos- 
mic star formation history has its broad peak. As a con- 



sequence, it should be possible to detect SNe la up to high 
redshifts, 2~ 5 so as to discriminate among different cosmo- 
logical models. In particular, at z>2 the SNe la rate should 
be nearly constant at a level of about lO"** SN yr~^ Mpc~^ 
(see Figure 6). Such a rate can be reduced only if the ef- 
fects of metallicity evolution become important at z~l, and 
if the changing in metallicity has an important effect in the 
explosion rate as predicted bv lKobavashi et alJ |l998^. 

• In the models predicting either bimodal or wide DTDs, 
"prompt" type la and Core-Collapse (CC) SNe are charac- 
terized by similar delay times and both trace the cosmic star 
formation history. At high redshifts the SN la "tardy" com- 
ponent tends to disappear and therefore we predict that the 
rate ratio CC/Ia steadily increases with redshift, as shown 
in Figure 8, from a value of about 3 in the local universe 
to about 9 at z^4. On the contrary, the "single-population" 
model in Figure 1 predicts a much faster evolution of the 
CC/Ia ratio, which is expected to become larger than 10 
already at z~1.5. 

Measuring the SN rates and their CC/Ia ratio at high 
redshifts will be a very interesting task for the upcoming 
James Webb Space Telescope and giant ground-based 
telescopes and will permit to confirm or discard these 
predictions. 
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